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Abstract In this paper we study the statistical properties of convex hulls 
of N random points in a plane chosen according to a given distribution. The 
points may be chosen independently or they may be correlated. After a non- 
exhaustive survey of the somewhat sporadic literature and diverse methods 
used in the random convex hull problem, we present a unifying approach, 
based on the notion of support function of a closed curve and the associated 
Cauchy's formulae, that allows us to compute exactly the mean perimeter 
and the mean area enclosed by the convex polygon both in case of indepen- 
dent as well as correlated points. Our method demonstrates a beautiful link 
between the random convex hull problem and the subject of extreme value 
statistics. As an example of correlated points, we study here in detail the 
case when the points represent the vertices of n independent random walks. 
In the continuum time limit this reduces to n independent planar Brownian 
trajectories for which we compute exactly, for all n, the mean perimeter and 
the mean area of their global convex hull. Our results have relevant applica- 
tions in ecology in estimating the home range of a herd of animals. Some of 
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these results were announced recently in a short communication [Phys. Rev. 
Lett. 103, 140602 (2009)]. 

Keywords Convex hull • Brownian motion • Random Walks 
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1 Introduction 

Convex sets are defined by the property that the line segment joining any 
two points of the set is itself fully contained in the set. In the physical world, 
convex shapes are encountered in many instances, from convex elements that 
ensure acoustic diffusion in concert halls [133], to crystallography, where the 
so-called Wulff construction leads, in the most general case, to a convex 
polyhedron whose facets correspond to crystal planes minimizing the surface 
energy [52]. Also, recent work in neuroscience indicates that the human brain 
seems better able to distinguish between two distinct shapes when these are 
both convex [80], which is particularly interesting since convexity properties 
are widely used in computer-aided image processing, in particular for pattern 
recognition [5]. Such applications would be limited if they were restrained 
to intrinsically convex shapes; but it is not the case, since it is possible to 
"approximate", in some sense, a non-convex object by a convex one: pick, 
among all convex sets that can enclose a given object, the smallest one in 
terms of volume. This is called its convex hull, and comparing convex hulls 
can be a viable mean of comparing the shapes of complex patterns such as 
proteins and docking sites [117]. Convex hulls thus attract much interest, 
both for the algorithmic challenges set by their computation [17, 51, 53, 74, 
83,88, 121, 142, 152, 158] and for their applications [5, 117, 145, 159, 160]. 

Random convex hulls are the convex hulls of a set of N random points 
in a plane chosen according to some given distribution. The points may be 
chosen independently each from an identical distribution, e.g., from a uniform 
distribution over a disk. Alternatively, the points may actually be correlated, 
e.g., they may represent the vertices of a planar random walk of N steps. 
For each realization of the set of points, one can construct the associated 
convex hull. Evidently, the convex hull will change from one realization of 
points to another. Naturally, all geometric characteristics of the convex hull, 
such as its perimeter, area, the number of vertices etc. also become random 
variables, changing their values from one realization of points to another. 
The main problem that we are concerned here is to compute the statistics 
of such random variables. For example, given the distribution of the points, 
what is the distribution of the perimeter or the area of the associated convex 
hull? It turns out that the computation of even the first moment, e.g., the 
mean perimeter or the mean area of the convex hull is a nontrivial problem. 

This question has aroused much interest among mathematicians over the 
past 50 years or so, and has given rise to a substantial body of literature 
some of which will be surveyed in Section 2. The methods used in this body 
of work turn out to be diverse and sometimes specific to a given distribution 
of points. It is therefore important to find a unified approach that allows one 
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to compute the mean perimeter and area, both in the case of independent 
points as well as when they are correlated such as in Brownian motion. The 
main purpose of this paper is to present such an approach. This approach is 
built on the works of Takacs [151], Eddy [54] and El Bachir [58] and the main 
idea is to use the statistical properties of a single object called the 'support 
function' which allows us, using the formulae known as Cauchy, Cauchy- 
Crofton or Cauchy-Barbier formulae [7, 11, 37,43, 137, 154], to compute the 
mean perimeter and the mean area of random convex hulls. This unified 
approach allows us to reproduce the existing results obtained by other diverse 
approaches, and in addition also provides new exact results, in particular for 
the mean perimeter and the mean area of n independent planar Brownian 
motions (both for open and closed paths), a problem which has relevant 
applications in ecology in estimating the home range of a herd of animals. 
The latter results were recently anounced in a short Letter [124]. 

Our unified approach using Cauchy's formulae also establishes an im- 
portant link to the subject of extreme value statistics that deals with the 
study of the statistics of extremes in samples of random variables. In the 
random convex hull problem when the sample points are independent and 
identically distributed (i.i.d) in the plane, then the associated extreme value 
problem via Cauchy's formulae is the classical example of extreme value 
statistics (EVS) of i.i.d random variables that is well studied, has found a 
lot of applications ranging from climatology to oceanography and has a long 
history [40,69,76,77]. For example, in the physics of disordered systems the 
EVS of i.i.d variables plays an important role in the celebrated random energy 
model [20,49]. In contrast, when the points are distributed in a correlated 
fashion, as in the case where they represent the vertices of a random walk, 
our approach requires the study of the distribution of the maximum of a set 
of correlated variables, a subject of much current interest in a wide range of 
problems (for a brief review see [110]) such as in fluctuating interfaces [32,78, 
79, 79, 105, 106, 122, 125, 139], in logarithmically correlated Gaussian random 
energy models for glass transition [36,63,64], in the properties of ground state 
energy of directed polymers in a random media [46,50,85,96, 108, 147] and the 
associated computer science problems on binary search trees [14,89, 109] and 
the biological sequence matching problems [112], in evolutionary dynamics 
and interacting particle systems [15,91, 111, 134,148], in loop-erased random 
walks [4], in qucucing theory applications [86], in random jump processes 
and their applications [39,41,107], in branching random walks [23,27,116], 
in condensation processes [59], in the statistics of records [70,90,97,115] and 
excursions in nonequilibrium processes [65,71, 146], in the density of near- 
extreme events [135, 136], and also in various applications of the random 
matrix theory [19,47,48,92, 104, 114, 119, 140, 153, 155]. Here, our approach 
establishes yet another application of EVS, namely in the random planar 
convex hull problem. 

The paper is organised as follows. In Section 2, we provide a non-exhaustive 
survey of the literature on random convex hulls. In Section 3, we introduce 
the notion of the 'support function' and the associated Cauchy formulae for 
the perimeter and the area of any closed convex curve. This section also es- 
tablishes the explicit link to extreme value statistics. We show in Section 4 
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how to derive the exact mean perimeter and the mean area of the convex hull 
of N independent points using our approach. Section 5 is fully devoted to 
the case when the points represent the vertices of a Brownian motion, a case 
where the points are thus correlated. Finally we conclude in Section 6 with 
some open questions. Some of the details are relegated to the Appendices. 



2 A (non-exhaustive) review of results on random convex hulls 

In this section we briefly review a certain number of results on the the convex 
hull of randomly chosen points. This review is of course far from exhaustive 
and we choose only those results that are more relevant to this work. They are 
presented in a chronogical order and at the end of the section we summarize 
in a table the results that are particularly relevant to the present work. 



In his book on random processes and Brownian motion (published in 
1948) [101], P. Levy mentions, in a few paragraphs and mainly heuristi- 
cally the question of the convex hull of planar Brownian motion: "This 
contour [that of the convex hull of planar Brownian motion] consists, 
except for a null-measure set, in rectilinear parts." 

More than ten years later, in 1959, J. Gcffroy seems to be the first to 
publish results pertaining to the convex hull of a sample of random points 
drawn from a given distribution [66], specifically N points chosen in K 2 
according to a Gaussian normal distribution /. He shows that if one 
denotes 

— by <9Cjv the boundary of the convex hull of the sample, 

— by En the ellipsoid given by the equation Q f(x, y) = jj, 

— by An the distance between 8Cn and En, 

— and by An the largest possible radius of an open disk whose interior 
lies inside the convex hull but contains none of the sample points, 

then almost surely: 

A N and A N — > (1) 



In other words, the convex hull of the sample "tends" to the ellipsoid 

N ' 

1961 [67]. 



riven by f(x,y) — -h-. Geffroy generalizes this result to R d , d > 1, 



1 For instance, in the Gaussian case: 

f(x,y) = exp 



2 i 21 

x +y 



and the ellipsoid is simply the circle centered on the origin with radius y/1 log N. 
We shall see further on how one can obtain directly the asymptotic behaviour of 
the perimeter of the convex hull in the case of n points chosen at random in the 
plane according to a Gaussian normal distribution. It is given by: 



(L N ) ~ 27rv/21ogJV, 



in complete agreement with Geffroy's result. 
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In 1961 F. Spitzer and H. Widom study the perimeter Ljy of the convex 
hull of a random walk represented by sums of random complex numbers 

S = 0, • • • , St = Z\ + Z 2 H Zk, 1 < k < N, where the Z k are i.i.d. 

random variables. By combining an identity discovered by M. Kac with a 
formula due to A.-L. Cauchy (which is used here for the first time in the 
context of random convex hulls) , they derive an elegant formula for the 
expectation [150]: 

E(L N )=2J2 E( ^, (2) 

fc=i 

The asymptotic behaviour of E(Ln) is studied in two different cases: 
1. Writing Z k = X k + iY k and taking E(X k ) = E(Y k ) = 0, E(X%) = a 2 , 
E(Y k 2 ) = b 2 , and E{X k Y k ) = pab, one has: 



E(L N ) ~ 4cViV, (3) 

where c(a, 6, p) does not depend on N. 
2. Taking Z k = X k + i with E(X k ) = p and E((X k - p) 2 ) = a 2 , one has: 



E(L N ) ~ 2Ny/l + ^ + - r-rbgiV, (4) 

which expresses the excess of E(Ln) over its smallest possible value 
2^71+^2 _ 

In the same year, and still for very general random walks viewed as a sum 
of TV vectors in the complex plane So = 0, • • ■ , S k = Z\ + Z^ + ■ ■ ■ Z k , 1 < 
k < N, G. Baxter [13] establishes three formulae involving respectively 

— the number Fn of edges of the convex hull of the random walk, 

— the number of steps Z k from the walk which belong to the bound- 
ary of the convex hull, 

— the perimeter Ljv of the of the convex hull . 

In the latter case, Baxter's formula coincides with Eq. ([2]), but Baxter's 
derivation rests purely on combinatorial arguments and does not make 
use of Cauchy's formula. Instead, it relies on counting the number of 
permutations of the random walk's steps for which a given partial sum 
S k belongs to the boundary of the convex hull. 
His formulae are: 

N 1 

E(F N ) = 2 V -~21og7V (5) 
z — ' m 

rn—l 

E(K N ) = 2 (6) 

k=l 

In 1963, appears the first [131] of two seminal papers by A. Renyi and 
R. Sulanke dealing with the convex hull of N independent, identically 
distributed random points Pi (i = I..N) in dimension 2. Denoting by Fn 
the number of edges of the convex hull, they consider the following cases: 



1. the Pi's are distributed uniformly within a convex, r-sided polygon K: 

E(Fjv) = ^r(log iV + 7) +T(if ) + o(l) (8) 

where 7 is the Euler constant and T(K) is a constant depending on 
K only and which is maximal for regular r-sided polygons, 

2. the P^s are uniformly distributed within a convex set K whose bound- 
ary is smooth: 

E(F N ) ~ a(K)N%, (9) 

TV— >oo 

with a(K) a constant depending on K, 

3. the Pi's have a Gaussian normal distribution throughout the plane: 

E(F N ) ~ 2V2^Tog7V (10) 

Af^oo 

The following year, 1964, the second of Rcnyi and Sulanke's papers [132] 
extend these results, focusing on the asymptotic behaviour TV — > 00 of 
the perimeter Ljy and area An of the convex hull of N points Pi drawn 
uniformly and independently within a convex set K of perimeter L and 
area A: 

1. If if has a smooth boundary: 

E(L N ) = L - 0{N~i) (11) 
E(Aat) = A-0{N~i) (12) 

2. If if is a square of side a: 

E(L N ) = Aa-0{N~^) (13) 

E(A JV ) = a 2 -0(^) (14) 

In 1965, B. Efron [57], taking cue from Renyi and Sulanke, establish 
equivalent formulae in dimension 3, together with the average number 
of vertices (faces in dimension 3), the average perimeter and the average 
area of the convex hull of N points drawn independently from a Gaussian 
normal distribution in dimension 2 or 3, or from a uniform distribution 
inside a disk or sphere: 

1. For instance, for N > 3 points in the plane, drawn independently from 
a Gaussian normal distribution, writing <j>(x) = (27r)~2 cxp(— \x 2 ) 
and $(x) = X!^ 4>{y) dy: 

E(V N ) = J°° 4> N - 2 (p)<t> 2 (p) dp, (15) 

E(L N ) = 4n Q J°° $ N -\p)<f(p) dp, (16) 



- DC 



E(A N ) = 3ttQ |°° ^"'WfeO dp, (17) 

where Vn, Ln and Ajy stand respectively for the number of vertices, 
perimeter and area of the convex hull . 



7 



2. In dimension 3, one has: 



E(F N ) 



^( N ) r ^(p^ip) dp, 



/ J —oo 



(18) 



E(E N ) 




(19) 



E(V N ) 



-E(F N ) + 2, 



(20) 



E(L N ) 




/ J —oo 



(21) 



E(A N ) 



12n( N ) r <S> N -\p)<p\p) dp, 



\ ^ J J -oo 



(22) 



F N and E N standing respectively for the number of faces and the 
number of edges of the convex hull (Ln is thus the sum of the lengths 
of the edges, and An the sum of the surface areas of the faces, Vn 
denotes again the number of vertices). 
He also computes the average volume of the convex hull of N vectors 
drawn independently from a Gaussian normal distribution (with zero av- 
erage and unit variance) in a space of dimension d < N: 



(for N = d + 1, this expression needs to be multiplied by 2). 

— In 1965 still, H. Raynaud communicates in the Comptes Rendus de 
I'Academie des Sciences [126], his generalization to R d of the formulae 
by Renyi-Sulanke and Efron pertaining to the number of vertices of the 
convex hull, cither in the Gaussian normal case or in the uniform case. 
In the case of a Gaussian normal case, with zero mean and variance a/2, 
Raynaud computes the probability density of the convex hull of a sample 
and shows that, in the limit when the number of points N in the sam- 
ple becomes very large, this distribution converges to a uniform Poisson 
distribution on the sphere of radius \J a log N centered at the origin. 

— In 1970, H. Carnal [35] addresses the question of the convex hull of N 
random points in the plane, drawn from a distribution which he assumes 
only to be circularly symmetric. He gives expressions for the asymptotic 
behaviour of the average number of edges, average perimeter and average 
area of the convex hull. In particular, he shows that the average number 
of edges, for certain distributions, goes to a constant (namely 4) when N 
becomes very large. 

— H. Raynaud publishes in 1970 a second paper [127] on the convex hull 
of N independent points (in both the Gaussian normal case throughout 
the space or the uniform case within a sphere) in R d . He gives detailed 
accounts of the results he announced earlier [126]. He also gives expres- 
sions for the asymptotic behaviour of the number of faces Fffl (or edges 
if d = 2) of the convex hull, and shows in particular that in the standard 



E(Voljv) 




Gaussian normal case: 



E(F ^ )) ^c O ^^ l0giV)|( ^ 1) (24) 



Note that for d = 2, one recovers Renyi and Sulanke's formula Eq.( 10 1. 
For d = 3, one has E(F^ 3) ) - ^TrlogiV 

Ten years later, in 1980, W. Eddy [54] introduces the notion of support 
function into the field of random convex hulls. Considering, in the plane, 
N points Pi — (xi,yi) with a Gaussian normal distribution, he associates 
to each a random process defined by: 

Bi(0) = Xi cos + yi sin 0, 

9 varying from to 2ir. Note that 23^(0) is simply the projection of point 
Pi on the line of direction 9. He further defines the random process 

M(0)=sup{A(0)}, 

i 

whose law is shown to be given by that of the pointwise maximum of the 
N independent, identically distributed random processes 23,(0) (cf [24]). 
Eddy then shows that the point distribution of the stochastic process 
M(9) is given by Gumbel's law, and he hints (without going further) to 
the fact that certain functionals of M{6) give access to some geometrical 
properties of the convex hull of the sample: 

L N = M(9)d9 (25) 

jO 

for the perimeter; and: 

A N = \ [Ad* (6) - (M'(0)) 2 ] d9 (26) 
z Jo 

for the area, where M'(9) = These functionals are what we refer to 
as Cauchy's formulae. 

It is precisely the first of these formulae that L. Takacs [151] suggests one 
should use to compute the expected perimeter length of the convex hull 
of planar Brownian motion, in his 1980 solution to a problem set by G. 
Letac in the American Mathematical Monthly in 1978 [98]. Denoting by 
L t the perimeter of the convex hull of Brownian motion 23(r),0 <r<t 
, Takacs shows that: 

E(L t ) = Vffri. (27) 

The calculation is performed using the support function of the trajectory, 
in the same way as Eddy [54] hinted at for independent points. Planar 
Brownian motion 2?(r) is written as (x(r), y(r)), where x and y are stan- 
dard one-dimensional Brownian motions. One then defines 



z T (9) = x(t) cos 9 + y(r) sin 0. 
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This stochastic process z T (8) being nothing else but the projection of the 
planar motion on the line with direction 9, it is itself, for a fixed 9, an 
instance of standard Brownian motion. Hence, the M(9) that appears in 



Cauchy's formula (Eq. (25 1) and which is defined as: 



M(9) = max{z T (6»)}, 

0<T<t 

follows for a given 9 the same law as the maximum of a standard one- 
dimensional Brownian motion. In particular, it is independent of 9 and 
therefore one can write: 

E(L t ) = 2tt E(M(0)) , 

thus using the isotropy of the distribution of planar Brownian motion. 
Knowledge of the right-hand part of this equation then yields the desired 
result. 

In 1981, W. Eddy and J. Gale [55] extend further the work started by 
W. Eddy [54] They point out the link between extreme-value statistics 
applied to N 1-dimcnsional random variables and the distribution of the 
convex hull of multidimensional random variables. They consider sam- 
ple distributions with spherical symmetry and distinguish between three 
classes according to the shape of the tails: exponential, algebraic (power- 
law) or truncated (e.g. distributed inside a sphere). Eddy and Gale then 
compute the asymptotic distribution of the associated stochastic process 
(the support function M(0)) when the number N of points becomes very 
large. The three classes of initial sample distribution yield three types 
of distribution for the limit process, given by the Gumbcl, Frechct and 
Weibull laws, which are well known in the context of extreme- value statis- 
tics. Eddy and Gale also remark that the average number of vertices of the 
convex hull in the "Frechet" case (that is, for initial sample distributions 
with power-law tails) tends to a constant (as proved by Carnal [35]). 
Following another route, N. Jewell and J. Romano establish the following 
year, in 1982, a correspondance between the random convex hull problem 
and a coverage problem, namely the covering of the unit circle with arcs 
whose positions and lengths follow a bivariate law [84]. Thus, for arcs of 
length 7r: 

Prob (circle covered) = Prob (conv. hull contains origin) 

and more generally, for arcs of lengths other than 7r: 

Prob (circle covered) = Prob (conv. hull contains a given disk) 

In 1983, M. El Bachir, in his doctoral dissertation [58], studies the convex 
hull C'(t) of planar Brownian motion B(t). In particular, he gives a proof of 
P. Levy's assertion [101] that almost surely C(t) has a smooth boundary. 
El Bachir also shows that C(t) is a Markov process on the set of compact 
convex domains containing the origin O. Denoting by dC(t) the boundary 
of C(t), he establishes: 

1. Prob(S(i) € dC(t)) = Prob(0 € dC(t)) = 

2. {t : B(t) £ dC(t)} has a null Lebesgue measure . 



El Bachir then computes explicitly, from Cauchy's formulae, the expected 
perimeter length and surface area of the convex hull of planar Brownian 
motion. For the perimeter, he derives a general formula for motions with 
a drif t it, of which the special case it = enables one to retrieve Takacs' 
V87rf. For the area, he obtains: 

E(A t ) = Y . (28) 

Over the following decade, the study of the convex hull of a sample of 
independent, identically distributed random points has attracted much 
interest. C. Buchta [28] has obtained an exact formula giving the average 
area of the convex hull of ./V points drawn uniformly inside a convex 
domain K, the existing formulae being so far mainly asymptotic. A few 
years later, F. Affentranger [3] has extended Buchta's result to higher 
dimensions, via an induction relation. Many details and references can be 
found in the surveys of Buchta [29], R. Schneider [141], W. Weil and J. 
Wieacker [157]. 

Another active route is the one open by Eddy [54] and Gale [55], whose 
works have been extended by H. Brozius and de Haan [26] to non-rotationally- 
invariant distributions. Brozius et al [25] also study the convergence in 
law to Poisson point processes exhibited by the distributions of quanti- 
ties such as the number of vertices of the convex hull of independent, 
identically distributed random points. The works of Davis et al. [45] and 
Aldous et al. [6] also follow this type of approach. 

Cranston et al. [42] resume the study of the convex hull C(t) of pla- 
nar Brownian motion and in particular of the continuity of its boundary 
dC(t). They show that dC(t) is almost surely C 1 and mention work by 
Shimura [143,144] and K. Burdzy [30] showing that for all a e (§,tt), 
there exist random times r such that C(r) has corners of opening a. They 
also mention Le Gall [95] showing that the Hausdorff dimension of the set 
of times at which the Brownian motion visits a corner of C (t) of opening 
a is almost surely equal to 1 — Finally, they also point to P. Levy's 
paper [102] and S. N. Evas' [60] for details on the growth rate of C(t), as 
well as to K. Burdzy and J. San Martin's work [31] on the curvature of 
C(t) near the bottom- most point of the Brownian trajectory. 
In 1992, D. Khoshnevisan [87] elaborates upon Cranston et a/.'s work by 
establishing an inequality that allows one to transpose, in some sense, the 
scaling properties of Brownian motion to its convex hull. 
In 1993, two papers concerned with the convex hull of correlated random 
points, specifically the vertices of a random walk, are published. G. Letac 
[99] points out that Cauchy's formula enables one to write the perimeter 
Ln of the convex hull of any iV-step random walk in terms of its support 
function Mn(8) — max <i<Ar{a:i cos 6 + yi sin 8}: 

/■27T 

E(L N ) = / E(M N (0)) d0, (29) 
Jo 

which provides an alternative to Spitzer-Widom's or Baxter's methods to 
compute the perimeter of the convex hull of a random walk. 
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It is precisely to the Spitzer-Widom-Baxter's formula (Eq. ([2])) that T. 
Snyder and J. Steele [149] return, using again purely combinatorial argu- 
ments to obtain the following generalizations : 

Let -Fjv be the number of edges of the convex hull of a iV-step 
planar random walk and let et be the length of the z-th edge. If / 
is a real-valued function and if we set Gn = J2f=i f( e i)> them 

EM^f^P, (30) 

fe=i 

where Sk — Z\-\-Z<i-\-- ■ ■ is the position of the walk after k steps. 

— Taking f(x) = 1, one has Gn — Fn (the number of edges of 
the convex hull) and one retrieves Baxter's result 

N 1 

k=l 

— Taking f(x) = x, one has Gn — Ln and one retrieves Spitzer 
and Widom's result without using Cauchy's formula: 

k=l 

— Taking f(x) = x 2 , Gn is the sum of the squared edges lengths 

(2) 

denoted by L N , one obtains: 



E 



(L^)=2iV(a| + 4), 



a\ + cry being the variance of an individual step. 
Snyder and Steele establish two other important results: 

1. An upper bound for the variance E (L N ^j (not to be mistaken for 

E (^L^^j) of the perimeter of the convex hull of any TV-step planar 
random walk: 

E(L%)<^N(a x + a 2 Y ) (31) 

2. A large deviation inequality for the perimeter of the convex hull of an 
TV-step random walk: 

Prob(|iiv -E(Ljv)| > t) < 2e~l^ (32) 

In a paper published in 1996 [72], A. Goldman introduces a new point 
of view on the convex hull of planar Brownian bridge. He derives a set 
of new identities relating the spectral cmpiral function of a homogeneous 
Poisson process to certain functionals of the convex hull. More precisely: 
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Let D(R) be the open disk with radius R and Di (i = l..N R ) 
the polygonal convex domains associated to a Poisson random 
measure and contained in D{R). Let 

oo 

Mt) = ^cxp(-iA njl ) 

n=l 

be the spectral function of the domain Di (the A^^'s being the 
eigenvalues of the Laplacian for the Dirichlet problem on Di). Fi- 
nally, set: 

Then: 

almost surely has a finite limit <£(i) (called the em- 
pirical spectral function) when R goes to infinity, and: 



^ = i E ( C ^ L ) (33) 



where L stands for the perimeter of the convex hull of 
the unit-time planar Brownian bridge (a Brownian motion 
conditionned to return to its origin at time t = 1). 
Goldman also computes the first moment of L using Cauchy's formula 



E(L) = A /- (34) 
To obtain the second moment, 

Goldman brings it down to computing E(M(#)M(0)), the two-point cor- 
relation function of the support function of the Brownian bridge, 



sin r 

E(M(0)M(O)) = — 



9(2tt - 9) 

— — + cotanfl 

6(7T — 0) 



(36) 



Goldman obtains this last result from the probability that the Brownian 
bridge £> ,i lies entirely inside a wedge £ of opening angle [3: 

Prober € = ^f> 2 (^) ^ 2 ), (37) 

with v = jf, and, assuming that O lies inside the wedge £, with r the 
distance between O and the apex 5 of the wedge, and with a the angle 
between the line OS and the closest edge of the wedge. (/„ is the modified 
Bessel function of the first kind.) 
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— In a later paper [73] , Goldman exploits further the link between Poisso- 
nian mosaics and the convex hull of planar Brownian bridges. He first 
shows that one can replace Brownian bridges by simple Brownian mo- 
tion. He then recalls Kendall's conjecture on Crofton's cell (in a Poisson 
mosaic, this is the domain D a that contains the origin): when the area 
Vo of this cell is large, its "shape" would be "close" to that of a disk. 
Goldman shows in this paper a result supporting this claim (in terms of 
eigenvalues of the Laplacian for the Dirichlet problem) and, thanks to the 
links he has established, deduces that the convex hull of planar Brownian 
motion, when it is "small", has an "almost circular" shape. More pre- 
cisely: if C is the convex hull of a unit-time, planar Brownian motion W, 
if D{r) is the disk centered at the origin with radius r e (0, oo) and if 
M = sup{||W(s)||,0 < s < 1}, then, for all e E (0,1): 



limsupProb[£>((l - e)a) C C C D(a)\M = a] = 1 (38) 



— In parallel, refined studies of the asymptotic distributions and limit laws 
of the number of vertices, perimeter or area of the convex hull of inde- 
pendent points drawn uniformly inside a convex domain K continue, in 
particular with the work of P. Groeneboom [75] on the number of ver- 
tices (augmented by S. Finch and I. Hueter's result [62]), Hsing [82] on 
the area when K is a disk, Cabo and Groeneboom [33] for the area too 
but when K is polygonal, Braker and Hsing [22] for the joint law of the 
perimeter and area, and the more recent works of Vu [156], Calka and 
Schreiber [34], Reitzner [128-130] and Barany et al. [8-10]. 

— In 2009, P. Biane and G. Letac [18] return to the convex hull of planar 
Brownian motion, focusing on the global convex hull of several copies of 
the same trajectory (the copies obtained via rotations). They compute the 
expected perimeter length of this global convex hull for various settings. 



Thus we see that random convex hulls have aroused much interest over 
the past 50 years or so. The main results relevant to our present study are 
those giving explicit expressions (exact or asymptotic) for the average perime- 
ter and area of the convex hull of a random sample in the plane. We have 
attempted to group the corresponding references in the following table: 
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Existing results 


Perimeter (average) Area (average) 


Independent 
Points 


Renyi and Sulanke [132] Eq. (fTTl and (13)); 
Efron [57] (Eq. (161 et (l7j}); Carnal [35]; 
Buchta [28[; Affentranger [3] 


Random 
walk 

(1 walker) 


open 
paths 


Spitzer et Widom [150] 

(Eq. @); 

Baxter [13]; 

Snyder et Steele [149]; 

Letac [99] 


? 


closed 
paths 


? 


? 


Brownian 
motion 

(1 motion) 


open 
paths 


Takdcs [151] Eq. (27)) 


El Bachir [58] 
(Eq. (281) 


closed 
paths 


Goldmann [72] 
(Eq. (341) 


? 



Finally, we have developed a general method recently [124] that enabled us 
not only to fill in the empty cells in this table but also to treat a generalization 
which is particularly relevant physically, namely the geometric properties of 
the global convex hull of n > 1 independent planar Brownian paths, each 
of the same duration T. To the best of our knowledge, this topic had never 
been addressed before. Among the works mentioned above, those dealing with 
correlated points always consider a single random walk or a single Brownian 
motion, except for [18] where several copies of the same Brownian path are 
considered. 

Yet, the convex hull of several random paths appears quite naturally, both 
on the theoretical side and also in the context of ecology, as we shall see later 



(§ 5.1 1. Furthermore, the simplest case, that of n independent Brownian mo- 



tions, already exhibits interesting features since the geometry of the convex 
hull depends in a non trivial manner on n [124]. We showed that even though 
the Brownian walkers are independent, the global convex hull of the union 
of their trajectories depend on the multiplicity n of the walkers in a nontriv- 
ial way. In the large n limit, the convex hull tends to a circle with a radius 
~ Vhi n [124] which turns out to be identical to that of the set of distinct sites 
visited by n independent random walkers on a 2-dimensional lattice [2,93,94]. 
This general method will be developed in detail in the following sections. 
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3 Support function and Cauchy formulae: 
a general approach to random convex hulls 

In this section we discuss the notion of the 'support function'. Intuitively 
speaking, the support function in a certain direction 9 of a given two di- 
mensional object is the maximum spatial extent of the object along that 
direction. We will see that the knowledge of this function for all angles 9 can 
be fruitfully used to obtain the perimeter and the area of any closed convex 
curve (in particular for a convex polygon) by virtue of Cauchy's formulae. 



3.1 Support function of a closed convex curve and Cauchy's formulae 

Let C denote any closed and smooth convex curve in a plane. For example C 
may represent a circle or an ellipse. The curve C may be represented by the 
coordinates of the points on it {(X(s), Y(s))} parametrized by a continuous 
s. Associated with C, one can construct a support function in a natural 
geometric way. Consider any arbitrary direction from the origin O specified 
by the angle 8 with respect to the x axis. Bring a straight line from infinity 
perpendicularly along direction 9 and stop when it touches a point on the 
curve C. The support function M{6), associated with curve C, denotes the 
Euclidean (signed) distance of this perpendicular line from the origin when 
it stops, measuring the maximal extension of the curve C along the direction 
9. 

M{9) = max {X(s) cos 9 + Y(s) sin 9} . (39) 

The knowledge of M(9) enables one to compute the perimeter of C and 
also the area enclosed by C via Cauchy's formulae 

L= d8M{9) (40) 
Jo 

A =lJ^ dd ( M2 ( d )-( M/ (8)f)- (41) 

These formulae are straightforward to establish for polygonal curves, and 
taking the continuous limit in the polygonal approximation yields the result 
for smooth curves (a non-rigorous, but quick, 'proof is provided in appendix 

In the elementary example of a circle centered on the origin with radius r 



(Fig. 1(a) I, M{6) is constant and equal to r for all 9. Its derivative is zero and 



Cauchy formulae give the standard results. The second example is slightly 



less trivial (Fig. 1(b) I as M{9) is not constant but equal to r(l + sin 9). One 



of course recovers again the usual results 

,-2n 



L= d9 r(l+sin0) = 2nr 
Jo 

1 r 27T 

A=- d9 r 2 [(l +sin60 2 - cos 2 <9] = irr 2 . 

2 Jo 
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It is interesting to note that Cauchy's original motivations for deriving 
the formulae (40 1 and (41 ) actually came from a somewhat different context. 
He was interested in developing a method to compute the roots of certain 
algebraic equations as a convergent series and to compute an upper bound 
of the error made in stopping the series after a finite number of terms. It 
was in this connection that he proved a number of theorems concerning the 
length of the perimeter and the area enclosed by a closed convex curve in a 
plane. Anticipating the usefulness of his formulae in a variety of contexts and 
particularly in geometrical applications, he collected them in a self-contained 
"Memoire" published by the "Academie des Sciences" in 1850. It is worth 
pointing out that Cauchy's formulae are of purely geometric origin without 
any probabilistic content. The idea of using these formulae in probabilistic 
context was first used by Crofton [43] whose work can be considered as one of 
the starting points of the subject of integral geometry, developed by Blaschke 
and his school during the years 1935-1939. 



3.2 Support function of the convex hull of a discrete set of points in plane 

Let I = {(xi,yi), i = 1,2,..., N} denote a set of TV points in a plane with 
coordinates (xi,yi). Let C denote the convex hull of /, i.e., the minimal 
convex polygon enclosing this set. This convex hull C is a closed, smooth 



convex curve and hence we can apply Cauchy's formulae in Eqs. (40 1 and 



(411 to compute its perimeter and area. However, to apply these formulae 



we first need to know the support function M{9) associated with the convex 



hull C, as given by Eq. (39). This requires a knowledge of the coordinates 
(A(s),y(s)) of a point, parametrized by s, on the convex polygon C. A 
crucial point is that one can compute this support function associated with 
the convex hull C of I just from the knowledge of the coordinates (xk,yk) 
of the set / itself and without requiring first to compute the coordinates 
(X(s),Y(s)) of the points on the convex hull C. Indeed the support function 



17 



M(6) associated with C can be written as 

M (9) = max {X(s) cos 8 + Y(s) sin 9} = max {x l cos 6» + ^ sin 6>} . (42) 



This simply follows from the fact that M(9), the maximal extent of the convex 
polygon C along 9, is also the maximum of the projections of all points of the 
set I along that direction 9. Thus, the knowledge of the coordinates (xk,yk) 
of any set / is enough to determine the support function M(9) of the convex 
hull C associated with / by Eq. (42 1. 




Fig. 2 Support function M{9) and its derivative M'(9) of the convex hull (dotted 
green lines) assocaited with a set of 7 points. 



We also note that, by definition of M(9) in Eq. (42), for any fixed 9 there 
will be a point (xk* , yk*) in the set such that: 



M (9) = x^ cos 9 + sin 9 



Taking derivative of ( 43 1 with respect to 9 gives 
M'{9) = -Xk* sin 9 + y k , cos 9 



(43) 



(44) 



In other words, M'(9) is the distance between the point of the set giving the 
maximal projection M{9) and the straight line with direction 9, as illustrated 
in Figure [2| 



3.3 A Simple illustration of the support function M(9) of a triangle 



To get familar with the support function M(9) associated with a convex hull, 
let us consider a simple example of three points in a plane. The associated 
convex hull is evidently a triangle (Fig . [3| w hose support function M(ff) and 
derivative M'(9) are drawn in Figures 4(a) and |4(b)| 

M(9) is of course 27r-periodic. A notable feature of its graph is the pres- 
ence of angular points, corresponding to discontinuities of the derivative of 
M{9). So, M{9) appears piecewise smooth, with a derivative exhibiting a 
finite number of jump discontinuities. This finite number is, in the special 
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Fig. 3 Triangle ABC (in red), with the line of direction 8 (in blue) and the line 
through O perpendicular to AB 



case considered here, equal to three, the number of points in the set whose 
support function is M(9). This is not by chance and the coincidence can be 



understood by returning to equations (143]) and ( 44 1 : within a given range 
of 9, one of the vertices of the triangle~A.BC will be giving the maximal 
projection on direction 9 and will thus determine the value of M(9), say: 

M(9) = x A cos 9 + y A sin9, (45) 

Then, within the same range of 9: 

M'{9) = -XA sin + y A cos 9. (46) 



We can specify the range of angles 9 on which equations (45 and 46 1 are 
valid. Let us indeed note that by the definition of M(9), when 9 corresponds 
to the perpendicular through the origin O to the line segment AB 7 A and 
B have the same projection on direction 9 (Fig. ~~). In the #~-limit, that 
is for angles approaching 9 from below, the value of M will be given by the 
projection of A, and the value of \M'(9) | by the length of the line segment AH 
(H being the foot of the perpendicular to AB through O). In the # + -limit, 
for angles slightly larger than 9, the value of M(9) will still be the common 
projection of A and B but |M'(0)| will be given by the length BH . Whence, 
as 9 passes on the perpendicular to AB through O, the support function M 
will be continuous, while its derivative will have a jump discontinuity. 

Indeed, looking at figure - ?} and starting from 9 = 0, we see that point A 
gives the maximal projection on direction 9, and the length of this projection 
decreases as 9 increases, until the direction given by 9 coincides with line 
(OH), which is the perpendicular to [AB] through the origin. At this point, as 
we have just noticed, M(6) has an angular point: B will then give the maximal 
projection, whose value will increase until 9 corresponds to line (OB) where 
M(&) attains a local maximum before decreasing until 9 coincides with the 
perpendicular to [BC] through O, where M(9) has a second angular point; 
and so on0 



2 In the specific example chosen here, all 3 vertices of the triangle are "visible" 
through a local maximum of M(0). However, this is not always the case. This 
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l 

0.8 
0.6 

M(6) 

0.4 
0.2 

-0.2 

(a) 



0.5 

M'(6) 



-0.5 



-1 

(b) 

Fig. 4 (a) Support function M(6) of triangle ABC and (b) its derivative M'(0) 



3.4 Cauchy formulae applied to a random sample 

Let us now examine how the Cauchy formulae can be applied to determine 
the mean perimeter and the mean area of a convex hull associated with a set 
of N points with coordinates (xi, j/j) in a plane chosen from some underlying 
probability distribution. The points may be independent or correlated. 



is easily seen by considering a configuration in which point H (the foot of the 
perpendicular to (AB) through O), while being by definition on the line (AB), is not 
on the line segment [AB]: for example, if// is beyond A, A will go "unnoticed". Yet, 
the coincidence of direction 9 with line {OH) will always result in a discontinuity 
of M'(9) (although without change in sign), which corresponds to an angular point 
for M{6). Thus the angular points of M(ff) count the number of sides (and, in 
dimension 2, of vertices) of the convex hull. As for the local maxima of M(6), they 
only count the number of vertices E of the convex hull that are such that the 
maximal projection on line (OE) is given by E itself — one could call such vertices 
"extremal" or "self-extremal" vertices. 
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For each i and fixed 9, let us define 

Zi (9) = Xi cos 9 + yi sin 9 (47) 
hi{6) = — a;, sin 6> + cos 0. (48) 

Zj is simply the projection of the i-th point in the sample on direction 9 and 
hi its projection on the direction perpendicular to 9. By definition (Eq. (Eg}): 



M(9) = miK{« i (fl)} = zt.(fl) (49) 

i 

for a certain index k* . 
One then has: 

M'{6)=h k *{d) (50) 



When the points (a;j, j/j) are random variables, so is the index A:* and sub- 
sequently both M(9) and M'(9) are also random variables. Taking averages 
in Cauchy's formulae ((|40|) and ((41 1) we get the mean perimeter and the 
mean area 

(L)= d9 (M(9)) (51) 
Jo 

<a> = \ C de ( {m2{6)) " ((M ' (e))2) ) (52) 

where (•) indicates an average over all realisations of the points, and we 
assume that this operation commutes with the integration over 9. 

In the most general setting, let us also define 

— He be the probability density function of the maximum of the Zi{9), i.e., 
of the random varibale Zk* (9) 

— pe be the probability density function of the index k* for which Zi{9) 
becomes the maximum 

— and Uifi be the probability density function of the random variable hi{9) 
for a fixed i and 9, 

then: 



(M(9)) - / z n e (z) dz (53) 

J — OO 

/oo 
z 2 (i e (z) dz (54) 
-oo 

((M'(9)) 2 ) = f f h? p e (k) a k , e (h) dk dh (55) 

J — OO 

p e (k) (hl(9)) dk (56) 
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With this formulation, it appears explicitly that random convex hulls are 
directly linked with extreme- value statistics, the study of extremes in samples 
of random variables. Indeed, when / is finite and the N points labelled by 
i £ I are chosen independently and are identically distributed, for instance in 
K 2 , then pg is the distribution of the maximum of N real-valued i.i.d random 
variables (namely the Zj(#)'s) — a classical example of EVS [40,69,76,77]. 
One can then use directly the results of the standard EVS of i.i.d random 
variables. On the other hand, when the points are correlated, we need to study 
the distribution of the maximum of a set of correlated random variables-a 
subject of much current interest as mentioned in the introduction. Here we 
need to go beyond i.i.d variables and take into acount the strong correlations 
between the random variables that changes the distribution of their maximum 
in a nontrivial way. 

If the sample points are the vertices of an TV-step 2-dimensional random 
walk, then the Zi(9)'s can be seen, for a fixed 9 as the vertices of an TV-step 
1-dimcnsional random walk, and pg is the distribution of the maximum of 
such a walk. Note that in this case, pg is the distribution of the step at which 
the 1-dimcnsional random walk Zi{&) attains its maximum [38,44,61,120]. 

One can also consider cases when / is not a discrete, finite set: e.g. the 
random set might be the trajectory B(t) = (x(t), y(r)) of a planar Brownian 
motion at times r € I = [0,T]. In such a case, both z T {9) and h T (0) are 
instances of 1-dimensional Brownian motion, and pg is the distribution of the 
maximum of 1-dimcnsional Brownian motion in [0,T], pg is the distribution 
of the time at which 1-dimensional Brownian motion attains its maximum 
in [0, T] (given by Levy's arcsine law [100]), and a T j the propagator of 1- 
dimensional Brownian motion between and r (i.e. the distribution of the 
position of a linear Brownian motion after a time r) . 

In the following section, we use this approach to compute the mean 
perimeter and the mean area of the convex hull of a set of N indepedently 
chosen points in a plane. In Section 5, we will examine how the same ap- 
proach can be adapted to compute the mean perimeter and the mean area 
of the convex hull of n independent planar Brownian paths each of the same 
duration T. 



4 Independent Points 

4.1 General case 

Let us consider here a sample of N points drawn independently from a bi- 
variate distribution: 

Prob (xj £ [x, x + dx],yi £ [y,y + dy}) = p(x, y) dx dy, 



Following the route explained in the previous section (Eqs. (47), (|48|), 
we let: 

Zi(9) = x, cos 9 + yi sin 9, 



22 



and: 



hi(e) 



-Xi sin 9 + yi cos ( 



4.2 Isotropic cases 

Let (xi, j/i), (x2, 2/2)7 • • • , (%N, Un) be TV points in the plane, each drawn in- 
dependently from a bivariate distribution p{x, y) that is invariant under rota- 
tion, i.e., p(x,y) — G(^f x 2 + y 2 ). In such an isotropic case, the distribution 
of the support function M(8) does not depend on 9 and it is thus sufficient to 
set = and hence Mx = M(0). The random variables Zi(0) and hi(0) are 
just, respectively, the abscissa x^ and ordinate yi of the points. Combining 
(51 1 and (53), we can then write the average perimeter of the convex hull 



(L N ) — 27r(max{x i }) = 2n(M N } 

i 

It is useful to first define the cumulative distribution 

F N (M) = Prob[Mjv < M]. 
For independent variables it follows that 

fM 



F N (M) 



Px{x) dx 



N 



(57) 
(58) 

(59) 



where px(x) = p(x,y) dy is the marginal of the first variable X. Thus, 
in a general isotropic case 



(M N ) = / M F' N (M) dM 

J — OO 

/OO 
M p x (M) 
-00 



2tt N 



Px(x) dx 

> 

M p x {M) Fjv-i(M) dM 



N-l 



dM 



(60) 



For the average area in the isotropic case, we can write it as (Eqs. (T52J) , 
(54} , @): 

(A^) = 7r(Af^)-7r(^.) ) (61) 

where y^* is the ordinate of the point (xk* , J/fe») with the largest abscissa, i.e. 
satisfying: 

x^ = max{xi} = -Mat. 



We can easily express the second moment of Mx that appears in (61 ): 

M 2 F' N (M) dM (62) 

= .V / M 2 p x (M) F N -i(M) dM. (63) 



(Mf 



N 1 
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To compute the second term in (611, that is, the second moment of the 
ordinate of the point with largest abscissa, we first compute the probability 
density function p of this point, which is defined by: 

Prob {(Xk*,Vk*) S [{x, y),(x + dx, y + dy)}} = p(x, y) dx dy (64) 
(Note that F N (M) = J p(M,y) dy.) 

It is not difficult to see that p(x k * , y k * ) can be expressed as the probability 
density that one of the N points has coordinates (x k *,y k *) and the N — 1 
other points have abscissas less than x*: 



p(xk> ,Vk*) = N p(x k . , y k * ) 



Px{x) dx 



Then: 



(k2. 



N 



y\, p(x k *,y k *) F N -\(x k *) dx k * dy k . 



(65) 



(66) 



It now suffices to insert (63 1 and (66 1 in (61) to obtain a general expression 
for the average area of the convex hull of N points drawn independently from 
an isotropic bivariate distribution p with marginal px- 



(A N ) = Ntt 



u 2 p x {u) F N -x(u) du 



N 7T 



v 2 p(u,v) Ftv-i(u) du dv (67) 



The equations (60 1 and (67) are the main results of this subsection. They 
provide the exact mean perimeter and the mean area of the convex hull of 
N independent points in a plane each drawn from an arbitrary isotropic dis- 
tribution. As an example, let us consider the case of a Gaussian distribution 
where the general expressions can be further simplified. Let 

P(x, y) 



2vr 



e 2 



We then have: 



and: 



Px{x) 



/2tt 



exp 



Px(x') dx' 



x') dx' = $(x) 



(68) 
(69) 
(70) 



Inserting these into equations (60) and (61 1, and performing suitable in- 
tegrations by parts, we obtain: 



L N ) = 4 7T 



{A N ) = 3 tt 



N-2 



(x) (j) 2 (x) dx 



N-3 



(x) 4> 3 (x) dx 



(71) 
(72) 



which coincide with the expressions derived by Efron [57] using a rather 
different method. 
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4.3 Asymptotic behaviour of the average perimeter and area 



To derive how the mean perimeter and the mean area behave for large N, 
we need to investigate the asymptotic large N behavior of the two exact 
expressions in Eqs. (60 1 and (61 1. For the mean perimeter, since it is exactly 
identical to the maximum Mjy (upto a factor 2ir) of N independent variables 
each distributed via the marginal px(x), we can use the standard analysis 
used in EVS, which is summarized below. For the mean area, on the other 
hand, we need to go further. We will give a specific example of this asymptotic 
analysis of the mean area later. 



Summary of standard extreme-value statistics 



Let zi, 22, • • • , Zn be independent, identically distributed random 
variables with probability density function p(z), and let = 
max re= i..jv {z K } be their maximum. Then 



N 

F N (M) = Prob (M N < M) 



M 

p(z) dz 

-oo 

In the limit when N becomes very large, the cumulative distribution 
function F^(M) exhibits one of the three following behaviours, ac- 
cording to the shape of the "tails" of the parent distribution p(z): 

1. When the random variable z has unbounded support and its dis- 
tribution p(z) has a faster than power law tail as z — > oo. We will 
loosely refer to this as "Exponential tails". Then, "Exponential 
tails lead to a Gumbel-type law" 

p(z) ~ -> F N (M) ~ e -e"<— 

z— +oo N— too 

2. "Power-law tails lead to a Frechet-type law" 

p(z) ~ A z-( Q+1 > -» F N (M) ~ e-- N M ~ a 

2^00 N^OO 

3. "Truncated tails (i.e. finite range a) lead to a Weibull-type law 
(with parameter a)" 

p(z) ~ A(a-z) a - 1 -» F N (M) ~ * (a-MT 

z—+a N^oc 

In all three cases, the typical value of the maximum Mjv increases as N 
increase^ the larger the number of points, the further the maximum is 
pushed. 



We will use these results for the general asymptotic behavior of the mean 
perimeter. However, before providing a summary of the asymptotic behavior 

3 As logiV in the first case, as a power of TV in the second, and nearing as an 
inverse power of N the radius a of the interval in the third case. 
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for a general isotropic distribution, it is perhaps useful to consider two special 
cases in detail, one for the mean perimeter and one for the mean area, that 
will illustrate how one can carry out this asymptotic analysis. For the mean 
perimeter, we choose the parent distribution from the Weibull-type case and 
for the mean area we choose the Frechet-type distribution. These choices 
are somewhat arbitrary, one could have equally chosen any other case for 
illustration. 

Example 1: average perimeter in the "Weibull-type" case 

Let us consider N points drawn independently inside a circle of radius 
a from a distribution with Weibull-type tails: 



p(x,y) 



A (a - y/x 2 + y 2 y-\ 



(73) 



Letting as before Fn denote the cumulative distribution function of 
the maximum Mjv of the x-coordinates, an integration by parts yields: 



(74) 



(M N ) = I x F' N (x) dx 
Fn(x) dx 



We focus on the second term of (|74|) and write: 
In-- 



Fn(x) dx. 



Then: 



IN 



exp 



Px{x') dx' 
N log [ I 



N 



dx 

Px{x') dx' 



dx. 



(75) 
(76) 



where as before we write px{x) = J p(x, y) dy. 
We now pick < e « 1 such that 

for (a — e) < x < a, 

p(x, y) ~ A (a - \J x 2 + y 2 ) 1 ^ 1 . 

The idea being that when N becomes large, some sample points will 
come closer and closer to the boundary (the circle of radius a) and con- 
sequently one can focus on the tails of the distribution. We therefore 
write: 



IN 



I 



(2) 
JV 



(77) 
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with: 



and 



r(2) _ 
1 N — 



exp 



exp 



N log 1 



Px(x') dx' 



N log 1 - / p x (x') dx' 



It is possible to show that iff decreases exponentially with N and is, 

(2) 

as expected on heuristic grounds, sublcading compared to I N which, 
as we are going to see, decreases as an inverse power in N. 

The sample distribution p(x,y) is rotationally invariant and bounded 
(x 2 + y 2 < a 2 ). Hence the marginal 



dx (78) 



dx (79) 



Px(x) 



V a 2 —x 2 



V a 2 —x 2 
V a 2 —t 2 



p(x,y) dy 



p(x,y) dy 



(80) 
(81) 



Now for x < a, we have: p(x, y) ~ A (a— x 2 + y 2 ) 1 1 . Consequently, 
setting y = xu and considering that (a — e) < x < a: 



\/a 2 — x 2 



Pjc(ss) ~ 2 

- 2Ax 



A {a~ ^x 2 + y 2 )~ f - 1 dy 



(a-x ^l + u 2 y- 1 du 



(a — x) 1 1 du 



2 A V2a (a-x) 7- 5 



(82) 

(83) 

(84) 
(85) 



r(2) _ 



exp 



exp 



We now proceed from equation ( 79 1 : 

N log^l- jf pjp^) dx' 
TV log ^1 - 

N log | 1 — 



exp 



exp - 



dx 

2 A V2a (a-x') 7_ 5 dx' 
4 A v/2a 



2 7 + 1 



(a-x^+i 



dx 



4 AN V2a 
2 7+1 



(a — x) 



7+1 



dx 



(86) 
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To progress further, we perform the following change of variable: 

A A N \[2a , +i 



2 7+ 1 



(a-x)- 



(87) 



In the large N limit in which we are working, this change of variable 
leads to: 



(2) 
N 



[4ANV2a] 1+2 -< j o 

2(2 7 + i)^ r(^ FI 



e~" [(27+ 1) du 



[AANy/^a] 1 



The combination of (|88| with (|74|) and (|57| yields the final result 
(Ln) 



2na 



47r(2 7 +l)^ 



[4ANV2a] M 



(89) 



To illustrate this asymptotic result for a concrete example, consider 
N points drawn independently and uniformly from a unit disk 



1 



0(1 -x 2 - y 2 ), 



where & is the Heaviside step function. 

In terms of our notations, this corresponds to: 

= 1, 

A-k 



7=1. 



We find: 



N 



N- 



2tt 



1 - 



7T3 



123 AT 3 



(90) 

(91) 
(92) 
(93) 

(94) 



in complete agreement with Renyi and Sulanke's result Eq. ( 11 1) [132]. 
Note that for large N, the mean perimeter of the convex hull ap- 
proaches 2n, i.e., the convex hull approaches the bounding circle of 
radius unity of the disk. But it approaches very slowly, the correction 
term decreases for large TV only as a power law ~ N~ 2 ^ 3 . Actually, 
for this example of uniform distribution over a unit disk, one can also 
obtain simple and explicit expressions for the mean perimeter and the 
mean area starting from our general expressions in Eqs. (60 1 and (61 ). 
Skipping details, we get 

Perimeter: 



(L N ) = 2tt 



1 - 



F N (M) dM 



(95) 
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Area: 



(A N ) = TT 1 - 8 [ M F N (M) dM 
3 J_i 



(96) 



where: 



F N (M) = - [arcsin(M) + M \J\-M 2 



N 



(97) 



One can also easily work out the asymptotic behavior of the mean 



area in this example using Eqs. (96 1 and (97 1 and we get 



(A N ) ~ 7r ll-2 r ® 2 '** ) (98) 



which, once again, agrees with Renyi and Sulanke's result (Eq. ( 12 1 ) 
[132]. Note also that the mean area of the convex hull approaches, for 
large N, to the area of the unit disk. Notice also that the exponent of 
N, which governs the speed of convergence is the same for the area as 
for the perimeter — only the prefactor of the power of N change^] 

At this point it is also worth recalling the results of Hilhorst et 
al. [81] regarding Sylvester's problem^] When N becomes large, the 
convex hull of the N points (conditioned to have all the N points 
to be its vertices) lies in an annulus of width ~ N~~s smaller than 
the N~3 found in our case. This can be understood qualitatively by 
noticing that requiring the N points to be on the convex hull will 
tend to increase the size of the hull and therefore push it closer to the 
boundary of the disk. 

Example 2: average area in the " Frechet-type" case 

Consider TV points drawn independently from an isotropic distri- 
bution with Frechet-type tails: 

p(x,y) - ^+2- (99) 

V^+y^oo (x 2 + y 2 ) 2 



Recalling equation (67 1, we start by its first term. Letting as before 
Fn denote the cumulative distribution function of the maximum Mjq 
of the x-coordinates, we have: 



(Ml) = / x 2 F' N (x) dx 



(100) 



4 Renyi and Sulankc [132] have shown that this is in fact true for every smooth- 

_ 2 

bounded support, and, moreover, with the same universal exponent: iV a . 

5 If N points are drawn from a uniform distribution in the unit disk, what is the 
probability pm that they be the vertices of a convex polygon — in other words that 
they be the vertices of their own convex hull? 
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With the same notation as previously: 



F N (x) 



/ px(x') dx 

J x 



N 



(101) 



where as before we write px{x) — J p(x,y) dy. 
We pick K ^> 1 such that for 

x>K, 

A' 

p(x,y)~ 



(x 2 + y 2 ) '' 2 2 



The idea being that when TV becomes large, sample points will dis- 
seminate further and further in the plane, and consequently one can 
focus on the tails of the distribution. We therefore write: 



with: 



and 



+/(2) 
1 N — 1 jv ' 1 N 



4 1} - / x 2 F' N (x) dx 

<J — CO 

iff = / x 1 F' N (x) dx 

JK 

r(i) 



(102) 
(103) 
(104) 

(2) 



It is easy to show that F N , as before, is sublcading compared to I y N ' . 
The sample distribution p{x 1 y) is rotationally invariant and so: 

/OO 
P{x,y) dy (105) 
-CO 

(106) 

Tpa". Consequently, setting 



= 2 / p(x,y) 
'o 



Now for x » 1, we have: p(x, w) ~ - 

(x 2 +y 2)^j 

y = ux and considering cases when x » 1 : 



Jo 

- 2Ar / 
Jo 



. /L^g d y 

o {x z + y z j 2 
1 



XT+ 2 (1+ U 2)- 

^+ir(i + i) 
c 

x^ 



du 



(107) 
(108) 

(109) 
(110) 
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where we have set C = Ay/n < . 

I 2 +1 j 

Now we can express Fjf(x) for large x and large n: 

-i N 



F N (x) 



1 



c 

7 



Whence, still for large x and large N: 



3-/7+1 
N 



dx' 



(111) 
(112) 



We now insert (1121 into (1041, setting u - 



nc . 

7 x~t ' 



r(2) 
'jV 



r l 



(113) 



Let us now examine the second term in equation (67). As before, we 
focus on the large- a; part of the integral, which will dominate. We 
rewrite it, so as to bring it down to the calculation that we have done 
in the previous paragraph: 



N 



OO POO 



K 



y 2 p(x,y) Fjv_i(a;) dx dy 



f°° y 2 p(x, y) dy 



K 
oo 



Px(x) 



F' N (x) x 



A 



K 



7-1 



C 

F' N (x) x 2 dx 



oo (l + U 2 )^ 



du dx 



(114) 



K 



This last integral is, up to the factor — the same as (1041; conse- 
quently, we obtain: 



(A 



N- 



1 - 



1 



7-1 



x 2 F' N {x) dx, 



K 



which, combined to (113 1 and simplified, yields: 

2 



(A 



N 



N- 



7 

7-1 



7 



r 2 l- 



(115) 



This coincides with the result of Carnal [35] . 
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Asymptotic results for general isotropic case: 

The large N asymptotic analysis for a general isotropic distribution, both 
for the mean perimeter and the mean area, can be done following the details 
presented in the above two examples. We just provide a summary here with- 
out repeating the details. 

Average Perimeter: 

— Exponential tails (p(x,y) ~ Ae~( x +y ) 1 when (x 2 + y 2 ) — ► oo) 

(L N ) ~ 2irlog 1/a N 

— Power-law tails (p(x,y) <~ A , 7+ iJ : 

(x 2 +y 2 y 2 ' 

where B(x, y) is the beta function. 

— Truncated tails (p(x,y) ~ A(a — \J x 2 + y 2 ) 7 ^ 1 ): 



with 



„ , (7+|)^r(r^) 

/(a, 7) = - 

(2AV2a)T+3^ 



Therefore, we do find as expected the distinction between the three differ- 
ent universality classes of extreme-value statistics. The sets of independent 
points drawn from distributions with exponential tails have, on average when 
TV becomes large, a convex hull whose perimeter increases more slowly (in 
powers of log N) than sets drawn from distribution with power-law tails (for 
which the growth of the perimeter is in powers of N) , which reveals the lesser 
probability of having points very far from the origin in exponential-tailed dis- 
tributions than in power-law tailed distributions. 

Average area: 

— Exponential tails: 

(A N ) ~ 7rlog° N 

— Power-law tails: 
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— Truncated tails : 

I A \ 2 A 8 /(a, 7) 
V 3 A^^+r 

with 

1-27 



/(a, 7) 



(2AV2a) 



We find for the area the same characteristics as for the perimeter as far as 
the relative growths of convex hulls are concerned, depending on the shape of 
the initial distribution of the points. It is particularly worth noting that, in 
the case of exponential-tailed distributions, the asymptotic behaviour of the 
average perimeter and area of the convex hull correspond to the geometrical 
quantities of a circle centered on the origin and with radius log" N (a being 
the characteristic exponent of the initial distribution's exponential tails) rj 



4.4 A non-isotropic case: points distributed uniformly in a square 

Let us now examine a non-isotropic case: computing the average perimeter 
of the convex hull of A" independent points distributed uniformly in a square 
of side a. The bivariate probability density of the sample can be written as: 

p(x,y)=\e( C ^-x 2 )0(~-y 2 ), (116) 
4 4 

where & is the Heaviside step function. 

We will consider, as described in the introductory part of this section, the 
projection of the sample on the line through the origin making an angle 9 with 
the x-axis. We write the random variable corresponding to the projection of 
a sample point z = x cos 9 + y sin 9. Its density will be given by : 

x 2 )Q(~-y 2 ) dxdy, (117) 



focus on < 9 < J and write 



)0( a --y 2 )dy. (118) 



This is to be compared with Geffroy's results [66], p. 4 



8(z — x cos 9 — y sin#)<9( — 



where 5 is the Dirac delta function. 

Using the symmetry of the square, we can 

z — y sin 9 



This enables us to simplify (1171: 
1 



?(*) = 



0( 



z — y sin 1 
cos 9 
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Enforcing the condition that any point of the sample lies inside the square 
and thus making the Heaviside functions non-zero, we have: 

(i) -§<!/<§ (H9) 

(h) -f < < f (120) 

••■ — ("!. £Z I^) < V < ^ (f , *±fc^) (121) 



and: 



- ^ (cos (9 + sin (9) < z < | (cos 6» + sin (9) (122) 



There will be 3 cases: 

1. For 

— - (cos + sin 9) < z < — (sin 9 — cos 9) , 

the y-coordinate will vary between — | and z+ J n C g Sg and therefore: 

z+ % (cos 6> + sin 6>) 
cos fc> sin 6* 

2. For 

- (sin 9 — cos 9) < z < - (cos — sin 0) , 

the y-coordinate will vary between — | and § and therefore: 

<?(*) = (124) 
a cos 

3. For 

- (cos — sin 0) < z < - (cos 9 + sin 0) 

the y-coordinate will vary between z ~ s ? n c ° sfl and | and therefore: 

, s §(cos0 + sin6>) - z 

9*=^ a . ' ■ 125 

cos 9 sin # 

To lighten the notation, let us write henceforth: 

a g = ^ (cos + sin 6») (126) 

b e = a 2 cos 9 sin 9 (127) 

Denoting as before by M^{9) the value of the support function of the 
sample at angle 9, that is, the value of the maximal projection on direction 
9, we have: 

(M N (9)) = f ° z Fg N (z) dz (128) 

J — ae 

= [zF'g. N (z)] a ° ae - r Fg, N (z)dz (129) 

J — ag 

= ae-I e , (130) 
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where: 



Fg, N {z) 
h 



-, N 



q(z') dz' 



— ag 
ag 



F 8< n(z) dz 



(131) 
(132) 



To compute Ig we make use of our knowledge of q(z) (Eqs. (1231, (1241, 



125} ) and we obtain: 



(M N (6)) =a g - 



cos 6 



sintftan^fl 



2 JV-i(27V + i) 2^(^+1) 



antf^-tan^ 1 . 



v /6 tan0 2 F 1 (^-A^;^ 1 ) , (133) 



3 tan 



2-Fi being a hypergeometric function. 

Using known facts about the asymptotic behaviour of hypergeometric series 
[1], (Mn(&)) can be seen to behave in the following way for large N: 



( M„ W )~a,-^ + „( ' ' 



This then yields the desired result: 



(L N ) = 8 / (M N (6)) 



4 a 1 



(134) 

(135) 
(136) 



This is the same as Renyi and Sulanke's [132], which they obtained from a 
different approach. Note that, as in the case of points distributed uniformly 
inside a disk, the average perimeter of the convex hull tends to that of the 
boundary of the support — here 4a, the perimeter of the square — when the 
number N of points becomes large. However, the convergence here is slower 
than for a support with a smooth boundary like the disk: versus . 

One can think that, physically and statistically, it is somehow "more difficult" 
for the points of the sample to reach inside the corners of the square, making 
the convergence of the convex hull towards the square all the more slower. 



5 Correlated Points: One or more Brownian Motions 

As mentioned before, one of the advantages of the support function approach 
that we use in this paper is its generality: it can be applied to samples with 
correlations as well as to samples of independent points. In this Section, we 
study, using this method, the convex hull of n planar Brownian paths, a topic 
that has so far been considered only in the n = 1 case [58,72,73,101,151]. 
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Beyond its interest from a theoretical point of view, the study of the 
convex hull of n planar Brownian paths can be motivated by a question of 
particular relevance to the conservation of animal species in their habitat, as 
we shall see before giving the details of results. 



5.1 Planar Brownian paths and home-range 

A question that ecologists often face, in particular in designing a conservation 
area to preserve a given animal population [118], is how to estimate the home- 
range of this animal population. Roughly speaking this means the following. 
In order to survive over a certain length of time, the animals need to search for 
food and hence explore a certain region of space. How much space one needs 
to assign for a group of say n animals? For instance, in the case of species 
having a nest to which they return, say, every night, the "length of time" is 
just the duration of a day. In ecology, the home range is simply defined as 
the territory explored by the herd during its daily search for food over a fixed 
length of time. Different methods are used to estimate this territory, based 
on the monitoring of the animals' positions [68, 159]. One of these consists 
in simply the minimum convex polygon enclosing all monitored positions, 
called the convex hull. While this may seem simple minded, it remains, under 
certain circumstances, the best way to proceed [21]. 

The monitored positions, for one animal, will appear as the vertices of 
a path whose statistical properties will depend on the type of motion the 
animal is performing. In particular, during phases of food searching known 
as foraging, the monitored positions can be described as the vertices of a 
random walk in the plane [12, 16,56]. For animals whose daily motion consists 
mainly in foraging, quantities of interest about their home range, such as its 
perimeter and area, can be estimated through the average perimeter and area 
of the convex hull of the corresponding random walk (Fig. p)). If the recorded 



Fig. 5 Convex hull of a 7-step random walk 

positions are numerous (which might result from a very fine and/or long 
monitoring) , the number of steps of the random walker becomes large and to 
a good approximation the trajectory of a discrete-time planar random walk 
(with finite variance of the step sizes) can be replaced by a continuous-time 
planar Brownian motion of a certain duration T . (fig. rol. 




o 
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Fig. 6 Convex hull of planar Brownian motion 



The home range of a single animal can thus be characterized by the mean 
perimeter and area of the convex hull of a planar Brownian motion of duration 
T starting at origin O. Both 'open' (where the endpoint of the path is free) 
and 'closed' paths (that are constrained to return to the origin in time T) 
are of interest. The latter corresponds, for instance, to an animal returning 
every night to its nest after spending the day foraging in the surroundings. 
As we have seen in our review of existing results, the average perimeter and 
area of the convex hull of an open Brownian path are known [58, 151], as is 
the average perimeter for a closed path [72] . It seems natural and logical to 
seek an extension of these results to an arbitrary number of paths (fig. [7]), 
both from a theoretical point of view and from an ecological one, since many 
animals live in herds. We show first how to use the support-function method 
for n = 1 planar Brownian paths and then for n > 1 . 



Fig. 7 Convex hull of 3 independent, closed Brownian paths, starting at the origin 
O. 



5.2 Convex hull of a planar Brownian path 

We consider here a planar Brownian path of duration T, starting from the 
origin O: 

B(r) = {x(r),y(r)) 
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with 

< t < T, 

x{t) and y(r) being standard 1-dimcnsional Brownian motions of duration 
T obeying the following Langevin equations: 

x( T ) = Vx(t) 

and 

where T) x (t) and r) y {r) are independent Gaussian white noises, with zero 
mean and rfe/to-correlation: 

(v.(r)v.(r')) = S(r-r'). 

Let us note incidentally that this implies: 

(x 2 (r))=r 

and 

iv 2 (r)) = r (137) 



Fix a direction 8. We use as before (Eqs. (47 1 and (48 1) the projection on 
direction 8: 

zq(t) — x(t) cos 8 + y(r) sin 8 

and 

hg(r) — —x(r) sin 8 + y(j) cos 8. 

Now, zg and h@ are two independent 1-dimcnsional Brownian motion (each 
of duration T), parametrized by 8. It thus appears that M(8) is simply the 
maximum of the 1-dimensional Brownian motion z$ (r) on the interval r G 
[0,7], i.e., 

M(8) = max IzaM]. 
re[o,T] L 

Furthermore, if we write r* the time at which this maximum is attained, 
then: 

M{9) = z b {t*) = x(t*)cos8 + y(T*)sm6. 
Deriving with respect to 8 gives: 

M'{8) = -x{t*) sm8 + y{r*) cos 8 = h g { T *). 

In words, if M(8) is the maximum of the first Brownian motion zg(r), M'(6) 
corresponds to the value of the second, independent motion hg{j) at the time 
-* when the first one attains its maximum, (c/. Fig. |8(a) and |8(b)] ). 



In particular, when 8 = 0, zo(t) = x(t) and ho(r) = y(r), and M(0) is 
then the maximum of x(t) on the interval r € [0, T] while M'(0) = y{T*) is 
the value of y at the time r* when x attains its maximum. 
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Recall that in isotropic cases, Cauchy's formulae (Eqs. (51 1 and (52l) 
simplify to: 

(L) = 2ir (M(0)) (138) 

(A)=7r(([Af(0)] 2 )-([M'(0)] 2 ))- (139) 

The planar motion that we are considering here is assumed to be isotropic 
and we will thus use this version of the formulae. 

The distribution of the maximum of a 1-dimcnsional Brownian motion 
x(t) on [0,T] is known, and given by the cumulative distribution function: 

F(M) = Prob[M(0) < M] = erf ( -= ) , (140) 



?2T 



with 

2 f z 2 
erf(z) = — = / e~ u du. 

V 71 " Jo 

The first two moments of this distribution are readily computed: 

(M(0))= " 

and 



([M(0)f)=T. 



Equation (138) then yields the average perimeter of the convex hull of a 



planar Brownian path, 



(L) = V8ttT. (141) 



It is slightly more complex to compute the average area enclosed by 
the convex hull as one then needs to compute ([Af'(0)] 2 ). Let us first re- 
call (Eq. fll37|) that for a given t* , 



E[y 2 (r*)]=r 
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(since y is a standard Brownian motion), the expectation being taken over 
all possible realisations of y at fixed r*. But r* is itself a random variable, 
since it is the time at which the first process, x, attains its maximum. One 
therefore also has to average over the probabability density of r* (which is 
given by Levy's celebrated arcsine law: pi(r*) = [r*(T — t*)]~ x / 2 /tt); this 
leads to: 

([M'(0)] 2 > = <*r*>=T/2; 



whence one obtains, via equation (1391, the exact expression for the average 



area enclosed by the convex hull of the motion: 



(A) = ~ (142) 



If we now consider a closed Brownian path in the plane, that is, one 
constrained to return to the origin after time T, the reasoning is completely 
similar, but for x(t) and y(r) which are now Brownian bridges of duration 
T: both start from the origin and are constrained to return to it at time T: 

x(0) = x(T) = 

y(0) = y(T) = 0. 

The distribution of the maximum of a Brownian bridge is also known, and 
its first two moments are given by: 

(M(0)) = 

and 

([M(0)] 2 > = \- 




Equation (138 1 gives us as before the average perimeter of the convex hull: 



(L) = ^ — . (143) 

To compute the average area enclosed by the convex hull y(r), let us first 
note that for a Brownian bridge, at a fixed time r*: 

..2/„*M_ T*(T-T*) 



T 



Let us then recall another well-known result: the probability density of t* is 
uniform. Thus, averaging on r*, with uniform distribution pi(r*) = l/T, we 
obtain: 

<Af'(0) 2 ) = (y 2 (r*)) = |. 

Finally, as before, equation ( |139[ ) leads us to an exact expession for the 
average area enclosed by the convex hull of a 2d Brownian bridge: 

(A) = ^- (144) 
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Results ( |141[ ) and ( |142[ ) had been computed by M. El Bachir [58 ] , usi ng 
the same approach as here, hinted at by L . Takacs [151]. Equation ( |143| is 
given by A. Goldman. The last result, (144), is, to the best of our knowledge, 
new, as are those in the next paragraph, regarding the convex hull of several 
planar Brownian motions. 



5.3 Convex hull of n planar Brownian paths 

As mentioned earlier, the method can then be generalized to n independent 
planar Brownian paths, open or closed. We now have two sets of n Brownian 
paths: Xj(t) and yj(r) (j = 1,2, ... ,n). All paths are independent of each 
other. Since isotropy holds, we can still use Eqs. ( 138| and (139}, except that 



Af (0) now denotes the global maximum of a set of n independent one dimen- 
sional Brownian paths (or bridges for closed paths) Xj(r) (j = 1,2, ... ,n), 
each of duration T, 

M(0) = max max [xi(t), x 2 (t), . . . , x„(t)] . (145) 
i<:?< n 0<r<T 

Let j* and r* denote the label of the path and the time at which this global 
maximum is achieved. Then, using argument similar to the n = 1 case, it is 
easy to see that M'(0) = ^(t*), i.e., the position of the j^-th y path at the 
time when the x paths achieve their global maximum. 

To compute the first two moments of M(0), we first compute the distri- 
bution P„[M(0), T] of the global maximum of n independent Brownian paths 
(or bridges) xj(t). This is a standard extreme value calculation. 



5.4 Open paths 

Consider first n open Brownian paths. It is easier to compute the cumulative 
probability, 

F n (M) = Prob[M(0) < M]. 
Since the Brownian paths are independent, it follows that 

F n (M) = [F(M)] n , 

where 

/ M 

F{M) = erf 



>2T 

for a single path mentioned before. 

Knowing this cumulative distribution F n (M(0)), the first two moments 
(M(0)) and ([Af(0)] 2 ) can be computed for all n. Using the result for (M(0)) 



in Eq. (138 1 gives us the mean perimeter, (ijv) = OL n \JT with 

/>oo 

a N = 4nV2^ / du u [erf (u)]™" 1 . (146) 
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The first few values are: 



ai = V87T = 5.013.., 
a 2 = 4-v/tt = 7.089.. 



a 3 = 24 



tan" 1 (1/V2) 



.333. 



(see Fig. 11 for a plot of a n vs. n). 

For large n, one can analyse the integral in Eq. ( 146 1 by the saddle point 
method giving: 

a„ ~ 27r-v/21ogn. (147) 
(Details of the analysis are given in appendix [C| 

This logarithmic dependence on n is thus a direct consequence of extreme 
value statistics [77] and the calculation of the mean perimeter of the convex 
hull of n paths is a nice application of the extreme value statistics. 



To compute the mean area, we need to calculate ([M'(0)] 2 ) in Eq. ( R39| . 
We proceed as in the n = 1 case. For a fixed label j and fixed time r: 

£[yf(T)] = T, 

which follows from the fact that yj(r) is simply a Brownian motion. Thus: 

E[yl(r*)}=r*. 

Next, we need t o ave rage over r* which is the time at which the global 
maximum in Eq. ( 145 1 happens. The probability density p n (i~*) of the time r* 
of the global maximum of n independent Brownian motions (each of duration 
T), to our knowledge, is not known in the probability literature. We were able 
to compute this exactly for all n (details are given in appendix [SJ. We find 
that 

where the scaling function f n (z) is given by 

2n 



/»(*) 



7T \J Z(l - z) Jo 



dxxe x [erf (x\fz)\ 



n-l 



(148) 



A plot of f n (z) for various values of n is given in Fig. (5.4|. 



It is easy to check that for n = 1 , it reproduces the arcsine law mentioned 
before. 

Averaging over r* drawn from this distribution, we can then compute 



([M'(0)] 5 



t*p„(t*) dr* 



Substituting this in Eq. ( 139 ) gives the exact mean area for all n, (A n ) = f3 n T 
with 



fi n = 4ny/ir 



du 



u [erf(w)]™ 1 ^ue " — g(u)j 



(149) 
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0-011 = 


1 


B-B n = 






5 




10 


S-i/D = 


5!) 



f„(t) 



1 



Fig. 9 Probability density f n (t) of the time t at which the global maximum of n 
Bro wnian motions, each of unit duration, is attained, as given by the formula in 
Eq. (1481). 



where 



9( u ) 



1 



e"" 2 /* dt 



2^ Jo ^/t{l-t) 
For example, the first few values are given by: 

0i = tt/2 



02 

03 



1.570.. 

7T = 3.141.. 

7T + 3 - V3 = 4.409. 



(Fig. 11 shows a plot of /3 n vs n). 

The large-n analysis (details in appendix [C]) gives: 



(150) 
(151) 
(152) 



0n 



2tt In n. 



(153) 



5.5 Closed paths 

For n closed Brownian planar paths one proceeds in a similar way. The dif- 
ferences are: 

— the cumulative distribution function of the maximum of a single 1-dimensional 
motion is not given by equation |140| but by: 

F(M) = 1 - e" 3 ^ 

— the propagator of the 1-dimensional motion obtained by projection on 
the x-axis is that of a Brownian bridge and so: 
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Fig. 10 Probability density g n (t) of the time t at which the global maximum of 
n Br owni an bridges, each of unit duration, is attained, as given by the formula in 
Eq. (pil. 



the probabilty density of the time t* at which the maximum of n 1- 
dimensional Brownian bridges occurs is given by (see details in appendix 
B): 

Pn{r*) = ~g n f ~ 

where the scaling function g n (z) is given by 
4n 



9n(z) 



u 2 e~* 



1 _ e -4«" z(l-z) 



du. 



(154) 



A plot of g n (z) for different n is given in Fig. (10 1. 

Following then the same route as for open paths, we find that the mean 
perimeter and area are given by: 



and 

where, for all n. 



(L n ) = a n (c)VT 
(A n ) = (i n {c)T 



3/2 n 



k=l 



n \ (_1)*+1 

k Jk 



ft, 1 1 7L 

,fc=l fc=2 



and 



/(*) 



Q (k - I)- 3 / 2 (fctan-^x/fc^I) - V^T) 



(155) 
(156) 
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The first few values are: 



and 



a x (c) = \A" 3 / 2 = 3-937. 

a 2 (c) = V^V^- 1/2) = 5.090.. 

a 3 (c) = V^(3/V2-3/2 + l/V6) = 5.732. 

/3i(c) = tt/3 = 1.047.. 

/3 2 (c) = 7r(4 + 37r)/24= 1.757.. 

/3 3 (c) = 2.250.. 



(see Fig. 11 for a plot of a n (c) and /3 n (c) vs. n). 



Large n analysis (details in section O shows that: 



x n (c) ~ 7rVzlrn 



and 



A»(c) 



■ In n . 



(157) 
(158) 
(159) 



(160) 



(161) 



smaller respectively by a factor 1 /2 and 1/4 than the corresponding results 
for open paths — as one's intuition might suggest, given that a closed path 
is enforced to return to the origin. 



5.6 Numerical simulations and discussion 



We illustrate our analytical results on the convex hull of planar Brownian 
motion with some elementary numerical simulations. These require first to 
generate Brownian paths, and then to compute numerically the convex hull of 
the paths. Here we have used a simple algorithm known as Graham scan [74j^] 

Let us recall the asymptotic behaviours of the exact formulae plotted in 
figure 11 which are given by equations (1471, (1531, (1601 et (161 1. 
For n open Brownian paths: 



(L n ) - 27iV2Tlnn 
(An) - 27rTlnn; 



and for n closed Brownian paths: 

(4 C) ) 
<^ c) > 



7iV2Tlni 

7T 

— Tlnn. 
2 



(162) 
(163) 



(164) 
(165) 



7 This algorithm is not the quickest one, but our aim was mainly illustrative. 
The question of convex-hull-finding algorithms is a classic one in computer science 
[51,142]. 
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Number of paths 



Fig. 11 Setting T = 1, t he an alytical results for average perimeter a n [Eq. (146 1, 
in red] and area f3„ [Eq. (149 1, in green ] of n open Brownian paths, and simil arly 
the average perimeter a„[c) [Eq. (155 1, in brown] and area /3„{c) [Eq. (1561, in 
purple] of n closed Brownian paths, plotted against n. The symbols denote results 
from numerical simulations (up to n = 10, with 10 3 realisations for each point). 



In both cases, the ratio between the average value of the area of the convex 
hull and the square of the average value of the perimeter takes, asymptoti- 
cally, the same value as in the case of a circle: 



(A n ) ^ l_ 
(L n ) 2 n-+oo 4-7T 



(166) 



Thus, hcuristically speaking, for large n, the convex hull of n planar Brownian 
paths approaches a circle, centered at the origin, whose radius is obtained 
dividing (L n ) by 2ir: 

' = ^ (167) 



R n = \/2Thn 



for n open paths; and: 



Tin n 



(168) 



for n closed paths. 

Note that for finite n the shape of the convex hull is far from being a 
circle. It is only in the n — > oo limit that it approaches a circle. Roughly 
speaking, a large number of trajectories smoothen their global convex hull 
into a circular shape. 



Let us also make another observation. In the limit of large n, the pref- 
actor for the average area of the convex hull of n open Brownian paths, 
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Fig. 12 Convex hull of 10 open Brownian paths. (Paths are independent and start 
from the origin.) 



namely log n, is identical not only to that of the average area of the con- 
vex hull of n independent points drawn each from a Gaussian distribution 

p(x, y) oc e s ; but also to that of the number of distinct sites visited by 
n independent random walkers on a lattice]^] Indeed the number of distinct 
sites visited by n independent walkers on a lattice {eg Z 2 ) has been stud- 
ied systematically by H. Larralde et al. [93,94] (see also [2]). For n lattice 
walks (n 3> 1) each of k step (steps are only allowed to neighbourings sites), 
Larralde et al. have identified three regimes according to the value of k. For 
the second of these regimes, the intermediate one, the system is in a sort of 
diffusive state so that the number of distinct sites visited (S n (k)} grows like 
the area of the disk of radius Vk, that is, proportionally to k, with a prefac- 
tor logn. In this regime, Acedo and Yuste [2] describe the explored territory 
as "a corona of dendritic nature [characterized by filaments created by the 
random walkers wandering in the outer regions] and an inner hyperspherical 



core [where there is much overlapping], (c/. fig. 13 1 



The detailed transposition between lattice-walk models and ours (n pla- 
nar Brownian motions of fixed duration T) requires much care but it is 
interesting to note that in the intermediate regime, the number of distinct 



We thank Hernan Larralde for drawing our attention to this point. 
9 In particular because of the transition to both continuous time (that is an 
infinite number of steps) and continuous space (a lattice constant that tends to 0). 
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Fig. 13 Figure from Acedo and Yuste [2]: "A snapshot of the set of sites visited by 
n = 1000 random walkers on the two-dimensional lattice. The visited sites are in 
white, the unvisited ones are in black and the internal gray points are the random 
walkers. The outer white circle is centered on the starting point of the random 
walkers and its radius is the maximum distance from that point reached by any 
walker at the time the snapshot was taken. The internal black circle is concentric 
with the former but its radius is the distance between the origin and the nearest 
unvisited site." 

sites visited by n independent lattice walkers grows like a circle with the 
same radius as that of the convex hull of n planar Brownian paths. 



6 Conclusion 

The method presented in this paper, based on the use of support functions 
and Cauchy formulae, allows one to treat, in a general way, random convex 
hull problems in a plane, whether the random points considered are inde- 
pendent or correlated. Our work makes an important link between the two- 
dimensional convex hull problem and the subject of extreme value statistics. 

We have shown here how this method can be implemented successfully in 
the case of independent points in the plane and when the points are correlated 
as in the case when the points represent the positions of a planar Brownian 
motion of a given fixed duration T. This method should be adaptable to 
treat other types of random paths, such as discrete-time random walks or 
anomalous diffusion processes such as Levy flights. 

In addition, we have shown how to suitably generalise this method to 
compute the mean perimeter and the mean area of the global convex hull of 
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n independent Brownian paths. Our work leads to several interesting open 
questions: 

— For example, can one go beyond the first moment and compute, for in- 
stance, the full distribution of the perimeter and the area of the convex 
hull of n independent Brownian paths? 

— For a single random walker of N steps, the average number of vertices of 
its convex hull is known from Baxter's work [13]: (-FW) = 2(14-1 /2 + 1/3+ 
... 4- 1/N) ~ 21og(7V) for large N. It remains an outstanding problem 
to generalise this result to the case of the convex hull of n independent 
random walkers each of step N. 

— Furthermore, we have only studied the convex hull of n independent Brow- 
nian paths. However, in many situations, the animals interact with each 
other leading to collective behavior such as flocking. It would thus be 
very interesting to study the effect of interaction between walkers on the 
statistics of their convex hull. 

— Finally, it would be interesting to study the statistics of the convex poly- 
tope associated with Brownian paths (one or more) in 3 dimensions. Let 
us remark that Cauchy's formula exist in higher dimensions and that it 
is possible to apply a method similar to the one dcvcloppcd here in order 
to compute the average surface areap°| 

Acknowledgements We wish to thank D. Dhar , H. Larralde and B. Teissier for 
useful discussions. 

A Proof of Cauchy's formulae 

We give here a quick "proof of Cauchy's formulatf^] 



We consider a polygonal curve and, without loss of generality, examine the 
integrals appearing in Cauchy's formulae on a portion of the curve corresponding 
to the configuration shown on figure |A| On the interval 9 G [—<j>i, fa], the value of 
the support function of the polygonaTcurve will be given by A. Writing R for the 
distance between the origin O and vertex A, we therefore have: 




(169) 




(170) 



M{6) = Rcosd 



(171) 



The first of Cauchy's formulae then gives: 




(172) 



details will be published elsewhere 

We thank Deepak Dhar for suggesting the idea of this demonstration. 
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Fig. 14 Portion of a polygonal curve (in red) near one of its vertices (A), with the 
lines through the origin perpendicular to the curve before A and after A (in green), 
together with the line of direction 6 (in blue) 



which is indeed the length of the curve (that is, H 1 A + AH2) between Hi and H2. 
As for the second of Cauchy's formulae, it gives: 




= — (sin<jf>2 cos 4>2 + sin <j>\ cos^i) 



which is indeed the area of the polygon OH 1 AH2- 

This proves the formulae for closed polygonal curves containing the origin. If 
the origin is outside, one can see that the signs of the various terms will lead to 
cancellations and the formulae will remain valid. Finally, taking the continuous 
limit yields the result for smooth curves. 



B Time at which the maximum of n 1-dimensional Brownian 
motion is attained 

Let us write: 

Prob(M„ = M, r* = r) = p„(r, M) dr dM, (173) 

where: 

— M n is the global maximum of the n Brownian motions, 

— r* is the time at which this global maximum is attained, 

— pn is the joint probability density function of M n and r*. 
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This probability can be written as the probability that one of the motions 
attains its maximum M at time r and the n — 1 others all have a maximum which 
is less than M: 

p n (r,M) = np 1 (r,M)[F(M)} n - 1 , (174) 

F(M) being, as before, the cumulative distribution function of the maximum of 
one standard Brownian motion on the interval [0,T]: 

F(M) = erf (^)' 



The joint probability density function pi(r, M) of the maximum M and the time 
r at which it happens can be computed using various techniques. The simplest of 
them is to use the Feynman-Kac path integral method, but suitably adapted with 
a cut-off [103,106]. This technique has recently been used [113,123] to compute 
exactly the joint distribution pi (r, M) of a single Brownian motion, but subject to 
a variety of constraints, such as for a Brownian excursion, a Brownian meander etc. 
The results are nontrivial [113] and have been recently verified using an alternative 
functional renormalization group approach [138]. For a single free Brownian motion 
(the case here), this method can be similarly used and it provides a simple and 
compact result 

Pi(r,M) = 3 e sr. 

7TT 2 \JT — T 



We then obtain the marginal distribution p n (r) by integrating out M. It has the 
scaling form p„(r) — A-/ n (~) where the scaling function f n {z) is given in Eq. 



(148 1 and is plotted, for various values of n, in Fig. (5.4 1 
For Brownian bridges, the reasoning is exactly t 



and pi differ: 



le same, but of course F(M) 



Pl (r,M) 



2T M 
Mr(T- r]j 



e 2t(T— r) 



F{M) = l-e 



(175) 



(176) 



which can both be derived using the technique we pointed to above [103, 105, 113, 
123]. 



We then obtain p n (j) = f^ffn (55) where 
, , 4n 



gn{z) = 



Jo 



2 - 

u e 



[ 1 _ e -*u' .(I-.)]' 



du. 



(177) 



Figure 10 shows a plot of g n (z) for different values of n. 



C Asymptotic behaviour 

C.l Open paths - average perimeter 

Letting M n be the maximum of n independent Brownian paths each of duration T, 
we write Prob(M n < M) = F n (M). This can be expressed in terms of the cumu- 
lative distribution function F\ (M ) of the maximum of a standard one-dimensional 
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Brownian motion: 



F„(M) = [F(M)\ 
2 



du 



u du 



(178) 



In the limit when n and M become very large, (1781 becomes 

2 



F n (M) ~ exp 



nln 1 



exp 



2 r 

/7T / J 



e " i In 

e~ u du 



/ft J M 

/2T 



Integrating by parts yields: 



du = 



x/2M 



O(l) 



Inserting this into equation (179 1, one obtains: 
F n (M) ~ e 



Here we write: 



5 = M - \/2Tln? 



(179) 



(180) 



(181) 
(182) 



assuming that <5 vanishes at large n, as we will be able to check a posteriori. It 
follows that: 



M = 2Tlnn 1 + 



\/2T In? 



2Tlnn 1 + 



2d 



\/2Thu 



Hence: 



M 2 - 2Tlnn ~ 2<5\/2Thm = 2\/2T Inn ( M - \/2T In 



Inserting this into ( 182 1 yields: 



F n (M) ~ e" 



"V 2 T " Var law) 



(183) 



Now, to compute the asymptotic behaviour of M n when n is large, let us start 
from: 



(M n ) = / MF' n (M) dM 
Jo 

~ / MF^(M) dM 

J A 



(184) 
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where A >> 1 will "disap pear " at a later stage. 
Combining this with Eq. (|183[), one obtains: 



(M n ) ~ / M 

J A 

Setting: 



one arrives at: 



^ e -^(M-vw^ )e - e -^W(--^) dM (185) 



2 log n 
T 



M - V2Tlnj 



T 



<M„) ~ / I V2rinn + y 



In the limit when n — > oo, this leads to 

(M„) r 



dy (186) 



V2Tlnne" ! 'e" e 
V2T In n 



(187) 



Thus, the average perimeter of the convex hull of n Brownian paths in the plane 
behaves for large n 



(L n ) ~ 27rV2Tlnn. 



(188) 



C.2 Open paths: average area 

We wish to compute the asymptotical behaviour (for n large) of the average area 
of the convex hull of n open Brownian p aths in the plane, all independent and of 
duration T. We apply Cauchy's formula ( 139 1 in the context of isotropic samples: 



(A n ) = 7T ((Mn) - ([M' n ] 2 )) 



(189) 



The first term of the r ight-hand side is easily computed (it suffices to substitute 
M 2 for M in Eq. |l85| ): 

(Mn)~2Tlan (190) 

To obt ain ( 153 1, it thus remains to show that (A n ) is dominated by (Mn)- 
Recall Eq. ((1491: 



(A n ) = 4nT s/n I duu [erf(u)] 



.») 



(191) 



where 



ff(«) = 



2>/7r Jo y/r(l-r) 



The integral in Eq. is dominated by the contribution from the large-w part. 

Therefore, we examine in this limit: 



1 



/ dr- . 

Jo ^t(1 - r) 



(192) 
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Setting t = 1 — y: 



f dr ; 1 = e~^~ = / rfy ' = e i-a (193) 
7o v/r(l-r) 7o Vy(T^) 



rfy ■ / e - 2 ( 1+ ^ (194) 
o V2/(l-2/J 

dy 1 e- y " 2 (195) 

o yy(i-y) 



We now write: 
This leads to: 



1 , 1 _^ _„2 2« d« 

(XT — - P 



_„= r 2^ 
7o « 2 



2e-" 2 P , 1 



dz — e 



— 2 

g — li /■£ 



u 

An integration by parts shows that the second term on the right-hand side is 
whence: 

1 1 _^ ^e-^ 



dr - =e~— ~ y. . (196) 

o x/r(l-r) 

Returning to Eq. ( |191[ ), where the integral is dominated by the contribution from 
the large u part, we have: 

{A n )=4nT^[ du u [erf(M)]"" 1 (ue'^ - ff (u)) (197) 



4nT^ J duu [erf(tt)]"- 1 Iwe"" - J (198) 

/•oo 2 

4nTV7r / duu 2 [erf(u)] n_1 e~ u (199) 
Jo 

n(M%) (200) 

27rTlnn. (201) 



C.3 Closed paths: average perimeter 

The calculation is analogous to that of subsection |C.1| indeed, the only differ- 
ence is that the cumulative distribution function F n of the maximum M of n 
one-dimensional Brownian bridges is given by: 



2M A 

F n {M) = [l-e-^] 



(202) 
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In the limit when n and M become very large, we then have: 
F n (M) ~ exp 



_ 2M^ 

n In 1 — e T 



2 l,.1_T , 



(203) 



We retrieve here the same equation as (1821, where T is replaced by We can 
thus deduce the result given in ( 1601: 



(L n ) ~ 7rV2T Inn. 



(204) 



C.4 Closed paths: average area 

We wish to compute the asymptotic behaviour (for n large) of the average area of 
the convex hull of n closed Brownian paths in the plan e, all independent and of 
duration T. T he analysis is similar to that of subsection |C.2| We apply Cauchy's 
formula (11391) in the context of isotropic samples: 



(A n ) = 7T ({M 2 n ) - ([M' n ] 2 )) 



(205) 



The first ter m on the right-hand side is, as before, easily computed from subsec- 
tions [03] and EU 

T 



(Ml) 



ln n 



(206) 



We now show that (A n ) is dominated by (M%). The equivalent of Eq. ( 191 1 for 
closed paths is: 



(A n ) = 



2nT 



du u 1 1 — e 



(l - e " ) (ue u - £?(«)) 



(207) 



where: 



g 4t(1-t) 

dr — =, 
>/t(1 - r) 



The integral in Eq. (207 1 is dominated by the contribution from the large-it part. 
Therefore, we examine g(u) in this limit. 



We now write: 



This leads to: 



i ri e ITIT^T 

3 = 8 y 



V 1 
4 Jo 



dr 



e ~ 4t(1-t) 

y^r^) 



(208) 
(209) 



r=--*. 
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Setting: 
then yields: 



v = 2uz, 



1 f 2 e !- 4 * 



1 f" dv e 



VI - 4z 2 4 J u n 

-it 2 ru 



1 e" 



4 u 
8u 



*e- (l+~) 



y/n erf(u) + O 



Returning to Eq. (2071, where the integral is dominated by the contribution 
from the large- u part, we therefore have: 



(An) = 







2nT 
2tcT 

m n. 



du u ( 1 — e 



2 \ ^ — 1 



(we u - ff(w)j 



j 2 /, -« 2 V 

dun 1-e I 



- 1 I -.=> ^e" u 
1 we 



8» 



du u ( 1 — e 



2 \ n— 1 2 



(210) 
(211) 

(212) 

(213) 
(214) 
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